Density-Matrix Approach to a Strongly Coupled 

Two- Component Bose-Einstein Condensate 

Andal Narayanan * 
Non-Accelerator Particle Physics Group 
Indian Institute of Astrophysics, Koramangala, Bangalore - 560 034, INDIA 

"a^^ 

ON . and 
ON 

Hema Ramachandran 

> _ 

O I Raman Research Institute, Sadashivnagar, Bangalore - 560 080, INDIA 

\o 

(N 
(N 



O 



The time evolution equations for average values of population and relative 



I phase of a strongly coupled two component Bose-Einstein condensate (BEC) 

I is derived analytically. The two components are two hyper-fine states, which 



are coupled by an external laser that drives fast Rabi oscillations between 
these states. Specifically, this derivation incorporates the two-mode model 
proposed in for the strongly coupled hyper-fine states |1, — 1) and |2, 1) 



' of ^^Rb. The fast Rabi cycles are averaged out and the rate equations so 



derived represent the slow dynamics of the system. These include the collapse 
and revival of Rabi oscillations and their dependence on detuning and trap 

X" 

5h ' displacement as reported in experiments of A procedure for stabilising 

vortices is also suggested. 



I. INTRODUCTION 

Observations of Bose-Einstein Condensation in trapped dilute alkali atoms have opened 
up both experimental and theoretical challenges to understand the properties of such sys- 
tems. The dynamical properties of a single condensate such as its collective excitations in 
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a trap due to a time dependent drive has been one such area of research 0], 0. The 
experimental reahsation of simultaneous creation and confinement of Bose- Einstein conden- 
sates (BECs) in several hyperfine states of a given species of atom [Q, 0, has led to 
investigations on dynamics of two or more overlapping condensates by coupling them with 
an externally applied laser field. In particular experimental realisation of binary mixtures of 
two hyperfine states namely |1, —1) and |2, 1) of^^Rb has established the following properties 
of this coupled condensed system 0, 

(a) These two states have magnetic moments which are same to the first order. How- 
ever, due to other small effects such as gravity, the nuclear magnetic moment and 
nonlinearity in the Zeeman shifts, the location of minima for the two states in the trap 
can be adjusted to be slightly different or exactly coincident. 

(b) Spontaneous inter-conversion from one state to the other is not seen due to the large 
difference in internal energies between these two states. The hyperfine energy is 6.8 
GHz. This makes the two condensates distinguishable. These can be selectively imaged 
by choice of an appropriate laser. 

(c) These condensate states possess a relative quantum phase which can be measured. 
This phase evolves with time the rate being proportional to the chemical potential 
difference between the two condensates. 

(d) An external laser drive couples these two systems and helps to coherently transfer 
population from one state to the other. 

The mixed condensates thus offer an ideal experimental apparatus to look for macroscopic 
realisations of dynamical effects like standard Josephson effects [Q. 

Theoretical calculations on such Josephson like oscillations in these coupled Boson 
Josephson Junctions (BJJ) ||10[, |[Tl| have shown several interesting dynamical effects. 
Recently however an experimental observation of an unexpected behaviour of these 
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coupled systems was reported. In the limit of sustained and large field strengths of the exter- 
nal coupling laser, that is when Q, the Rabi frequency, was five to ten times larger than the 



trap frequency in the vertical direction, along which the two condensates sit displaced 12 



the Rabi oscillations between the hyperfine states was found to collapse and revive. This 
occurred on a time scale which is large compared to the Rabi period. These slow varying 
modulations of the fast Rabi oscillations vanish at zero trap displacement. It was also seen 
to vanish when 6 = 0, where 6 = uj - ujd is the detuning of the external laser frequency (ud) 
from the transition frequency between the hyper-fine states (a;). It was shown subsequently 
in the same paper that this phenomenon was due to a weak coupling between the low lying 
motional states of the trap. In particular a simplified two-mode model was suggested. In 
this two-mode model the trap ground state and first excited dipole state were coupled and 
couplings to all higher motional states were neglected. Thus the paper [0, demonstrates 
the possibility of quantum state engineering of topological excitations, through the interplay 
between the internal and motional degrees of freedom of a BEC in a TOP trap. Numeri- 
cal simulations by solving the Gross-Pitaeskii (GP) equations for the coupled system were 
carried out in which reproduced the experimental features. 

In the present paper, we derive the essential experimental features analytically, using 
the density matrix approach. Equations for the fractional population (Z) in the hyperfine 
states and their relative phase {9) as a function of time have been obtained. Averaging over 
the Rabi period, these equations represent the slow dynamics of the system. Collapse and 
revivals of Rabi oscillations and their dependence on detuning and trap displacement are 
seen to match qualitatively with the experimental results described in A proposal for 
combining this strongly coupled regime to the weakly coupled Josephson regime is presented 
and its role in increasing the stability of vortices ill], [ITH], [Il6l is also speculated. 
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II. A DENSITY-MATRIX METHOD 



The total wave-function of the two- component condensate is denoted by '0(r, t). Initially 
this wave-function just represents the total population (Nt) in the ground state |1, —1). An 
external laser drives the transition from this state to the |2, 1) state coherently. So we can 
write, 

7/;(r,t) = V^i(r,t) + ip2ir,t) 

These two states are of the form 

l^i) = (ai(t)co(t) |0o) + a2{t)d^{t) |1) (1) 
1^2) = («2(t)co(t) |0o) + alit)diit) |2) (2) 

Here |1) and |2) refer to the hyperfine (internal) states and |0o) and refer to the motional 
(external) states and 

«i(0=cos(%^)-z(-^)sin(%^) 
«2(t) = -^(7^)sin(^) 



fie// = V52 + fi2 

co(t) = COS (^) - H-^) sm (^) 
di{t) = -H — ) sm (— -) 



< z >ij= / (piZ(j)jdz 



Aeoi = ei - Co 



Here 6 is the detuning, Aeoi is the energy difference between the two trap states, namely 
the ground state and the first excited dipole state. This energy difference is held fixed 
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through the derivation, while in actuahty they will vary with time, zq is the displacement 
between the two condensates and < z >ij is the dipole matrix element which couples the 
ground and excited states of the trap. In this derivation < z >ij is held fixed at < z >oi. 
The higher couplings are weak and hence neglected. 

Taking equations ( [l|) and ( ||) as starting points, we point out that the and \4'2) 
individually satisfy the following normalized coupled Gross-Pitaeskii (GP) equation in the 
Thomas-Fermi limit in an isotropic trap. 



dt i 



= -{[-2zo < z > +X,N2 + iVi - ^]Mt) + -Mt)} (4) 



Nt = Ni + N2 
a± - 



I muj 

a = ail ~ O22 

, On 
Ai — 

ai2 

\ _ ^22 
A2 — — 

ai2 

In writing the above set of coupled equations, time is in units of trap frequency uoz and the 
spatial variables are scaled by a^. an and are respectively the s wave scattering lengths 
for the two hyperfine species of the condensates and ai2 is the inter-species scattering length. 
The various energy terms are given with respect to the trap energy level huo^. Due to a near 
degeneracy of an and a22 scattering lengths in the case of ^'^ Rh the approximation that they 
are equal can be safely carried out. The system is characterized then by a single scattering 
length a and a single A = Ai ~ A2. In deriving these equations, the spatial dependence of the 
GP wave-functions are integrated out (Adiabatic approximation) with respect to the trap 
wave-functions, namely \(j)Q{z)) and |0i(2;)) and treated as constants. This assumes that the 
specific changes in the shape of the trap wave-function is not playing a major role in the 



time evolution of the system. This assumption is an approximation over what is actually 
experimentally seen [Q] since we are only interested in the dynamics of the fractional 
population of the hyperfine states. 

We proceed to derive the population fractions and the relative phase differences between 
the two hyperfine states by noting that 



|l) = yiVi(t)e^^iW (5) 



|2) = v/^e^^^W (6) 

e = Lp2- (8) 

(9) 

By taking appropriate inner products of and |^/'2) with |0o) and |0i) and substituting 
for l-i/^i) and \il)2) in the GP equations, the form given in equations ( ^ and ( 0), the following 
rate equations for average values of fractional population {Z) and relative phase {9) in the 
two hyperfine states can be derived. 

< ^ >2 \ 



01 



■ sin {Vtmt) Z (10) 



(3^<_z >2 



^) = 4(Aeoi^^^^f— sin2(noit)__^^_ _ 



cos2(^) + (^)2sin2r^^ 



(11) 



In deriving the above equations the orthonormal relations of trap wave-functions are assumed 
to be 



I 



) = ^ij 



Equations ( p!0| ) and ( |TT]) are obtained after averaging over the fast time period namely that 
of Q in the problem. So these equations do not explicitly contain Q. Analytical expression 
for (Z) can then be derived. 

{Z) = Zoexp [^-^^ [1 - cosfioit] j (12) 



In the limit of small detuning S {5 '^Q), this is of the form 



1 + 2 -2 sm (— t) 



(Z) = J 

Here Zq is the initial value of the population at time t = 0. 



(13) 



III. RESULTS AND DISCUSSIONS 

The equation ( [T^) has all the essential features which are reported by the experiments 
and subsequent numerical investigation in [|l|. 

(1) The (Z), remains a constant when 5 goes to zero or when zq goes to zero. That is, in 
the lab frame the fast Rabi oscillations remain unmodulated. 

(2) Equation ( |13[) is derived with the implicit assumption that the condensate has a well 
defined overall phase which can be measured relative to a reference. A decoupling of 
{0) and (Z) in the time averaged frame over the fast variable occurs as this phase 
(both the slow and fast varying part) averages to zero. 

(3) Though the mean- field term does not explicitly enter the expression ( |13|), we can 
see that the amplitude of modulation increases with decreasing Acqi - a result which 
is confirmed by numerical simulations in [|l| which predicts a decrease in Acqi for 
enhanced mean-field effects. 

(4) This form of equation ( |1^) does not give rise to the chaotic behaviour with high values 
of Zq reported in [Q]. 

Figure 1 gives a typical curve for the parameters given. 



In the paper [Q, a preparation of the vortex mode is presented in this very two- 
component system. In this contest it is tempting to think of the following scheme for 
stabilising such vortices. Starting from the strong coupling regime (f2 > u^), the strength 
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of the external laser field could be gradually decreased over time. At sometime then, when 
Q < uJz, the population in a particular state (which is a combination of motional and 
internal states) gets trapped in that state itself, due to Macroscopic Quantum State Trapping 
(MQST) effects [|1T|, applicable in this regime. More specifically, if such trapping should 
occur in the first excited motional state, which could be a vortex state, then there seems to 
be a tremendous improvement achieved in the stability of the vortex. A schematic of such 
a scenario is given in Figure 2. 

IV. CONCLUSIONS 

In this paper we have analytically derived an expression for the rate of change of fractional 
population of the hyperfine states |1, —1) and |2, 1) of ^''^Rb in the strong coupling regime 
using the density matrix approach. The derivation gives analytical results for population 
evolution after averaging out the fast dynamical variable namely the Rabi period in the 
problem. This derivation is based on a two-mode model for the trap states as proposed in 
The main result of our analytical approach is presented in Equation ( |T3|). This equation 
reproduces most of the essential features of the two-mode model presented in Also the 
possibility to increase the stability of vortex state by modulating Q is also discussed. 
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TABLE I. Parameter values used to plot Figure 1 are given here 



Nt 


10^ 




65 Hz 


ail 


1.0 021 


0_L 


1.3 //m 


(122 


1.0 021 


021 


5.5 nm 
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1.00 ^ 
0.90 
0.80 ^ 



0.50 
0.0 



1.0 

time in units of trap frequency 




time in units of trop frequency 

FIG. 1. This plot shows the change in (Z) as a function of detuning 5. The parameters are (a) 
5 = 0, (b) 5 = 27r * 50Hz, (c) (5 = 27r * lOOHz. The values of other parameters are given in Table 1. 
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Schematic illustrating creation of stable vortex state 
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Figure 2. 
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